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I. INTRODUCTION 



There have been a number of significant proposals for 
new optical clocks that are more accurate than the cur- 
rent standard; the Cs clock 0,113. Singly charged ions 
are some of the prominent candidates in this category 
due to the remarkable advances in modern ion trapping 
and laser cooling techniques fil 0. H, 0| . S-D transitions 

in Sr+ d H, Yb+ H> H S + EH' Ba+ B S and Ra+ 
[Tfj| can serve as clock transitions. It is necessary to es- 
timate the shifts of the energy levels for these cases due 
to stray electromagnetic fields in order to determine the 
accuracies of these potential optical clocks. A knowledge 
of hyperfine structure constants, electric quadrupole mo- 
ments, gyromagnetic constants and polarizabilities are 
essential quantities that can be used to estimate various 
possible shifts [Tj| . In our previous studies, we have al- 
ready calculated hyperfine structure constants and elec- 
tric quadrupole moments for the above ions mm. it 
is possible to find the gyromagnetic constants approx- 
imately using analytical approaches for different states 
(TH . There have been extensive studies of the hyperfine 
structure constants of the low-lying states in the above 
ions using the relativistic coupled-cluster (RCC) method 
[T5L \l§L [lj, [ll| . A few calculations of polarizabilities in 
these ions using the sum-over-states approach have also 
been reported recently [h| EH • 

Both Ba + and Ra + have also been proposed as suit- 
able candidates for atomic parity violation (APV) experi- 
ments [171 Il8j . Determination of polarizabilities depends 
on the electric dipole (El) matrix elements and excitation 
energies. On the otherhand, the determination of APV 
amplitudes also depends on El matrix elements and ex- 
citation energies. Therefore, studies of correlation effects 
in these properties involving quantities in these systems 
are also useful for the APV studies. In contrast to hy- 



perfine structure constants where the explicit behavior 
of electron correlation has been studied elaborately, the 
same cannot be done for polarizabilities using the sum- 
over-states approach. Also, the sum-over-states approach 
considers only a limited number of states, mainly from 
the single excited states and misses out contributions 
from continuum, double excited states, normalization of 
the wave functions etc. 

In this work, we have employed an ab initio method 
in the RCC framework to calculate dipole polarizabili- 
ties of Sr + , Ba + and Ra + . The roles of different types 
of electron correlation effects in determining these quan- 
tities are studied and comparisons between these results 
are given explicitly. Contributions arising from differ- 
ent types of excited states and normalization of the wave 
functions through the RCC method have also been eval- 
uated. Using these results, we then estimate the Stark 
and black-body radiation shifts in these systems which 
will be useful for the proposed optical clock experiments 
in the ions mentioned above. 



II. THEORY 

The polarizability of a given state \ J n M n ) can be ex- 
pressed by 



2 £ e , M (-i)^ ' J -"^-'W |n "ft 



where the subscript represents for the static values and 
angular momentum coefficients (Cj(J„)) for the scalar 
(with superscript 1) and tensor (with superscript 2) 
dipole polarizabilities are given as 
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respectively. Es are the energies of the corresponding 
atomic states. 

By defining a modified wave function due to the dipole 
operator, D, we can rewrite the above expression as 



4(J n ) = (J B ||A||4 1) > + <4 1) HA||.7n>, (2.4) 



where 



i^) = E \jJ Jm]mJn) 



E n — E r , 



(2.5) 



appears as a first order correction to the \J n ) state due 
to the dipole operator D. In the above expression, we 
also define an effective dipole operator as 



Di 



Ci(J n )(-i) J »- J »A 



(2.6) 



whose matrix element between the original and per- 
turbed wave functions will give the dipole polarizabilities. 
Here the intermediate states, |J m ), have parities opposite 
to that of \J n ) and they have to satisfy the usual trian- 
gular condition for the vector operator D. 

To avoid the sum-over-states approach in the determi- 
nation of the polarizabilities, we avoid the explicit form 
of \J n 1] ) that is given by Eq. (|2.5| . In stead, we obtain 
| j}i ) by solving the following equation 



(H-E n )\jU) = -D\J n ) 



(2.7) 



that is similar to the first order perturbative equation. 
Here H is the atomic Hamiltonian which in the present 
work is considered in the Dirac-Coulomb approximation 

h = J2 [«* ■ Pi + - i y 2 + Vnucin)] + —.(2-8) 

i i>j 

where c is the velocity of light, a and (3 are the Dirac 
matrices and V nuc (r) is the nuclear potential. 



For an atomic system with zero nuclear spin, the Stark 
shift to the second order in the presence of an electric field 
(quadratic Stark shift) for | J„, M n ) state is given by [l9| 



,(2.3) AW E (J n ,M n ;E) 



-\al{J n )E 2 - \al(J n ) 



[ZMl - Jjjn + 1) 
J n (2(/ n 1 ) 



■(3.E; - E 1 



(2.9) 



where E and E z are the magnitudes of the externally 
applied electric field in any arbitrary and z directions, 
respectively. 

For atomic systems with non-zero nuclear spin (J), the 
expression for hyperfine states are given by 



AW E (F n ,M Fn ;E) 



-\a\(F n )E 2 - \al(F n ) 



[3M| 



F n (F n + 1)] 



F n (2F n - 1) 



(3E 2 - E 2 



(2.10) 



where F n = I + J n and Mp n are the total spin due to 
nuclear spin / and atomic state spin J n and its azimuthal 
component, respectively. Since it is easier for us to deal 
with J n of the electronic states, therefore we express all 
the above quantities in terms of electronic coordinate. 
By using the following relations [ll[ 



(2.11) 



and 



a 2 {F n ) = (-1) 



I+J n +F nn 2 



a (J n 



(2J„ + l)(J„ + l)(2J„ + 3) 



J n (2,J n 1 ) 
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between the dipole polarizabilities of the electronic and 
hyperfine states, we obtain 



AW E (F n ,M Fn ;E) 
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Again, the blackbody-radiation (BBR) shift of a given state | J„, M n ) in the adiabatic expansion due to the ap- 
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plied isotropic electric field radiated at temperature T (in 
Kelvin (K)) can be assumed as [20| 



*BBR 



i(831.9V7m) 2 al(J n ). (2.14) 



III. METHOD OF CALCULATIONS 

The RCC method, which is equivalent to all order per- 
turbation theory, has been recently used to obtain precise 
results and account for the correlation effects accurately 
in single valence systems 0, El EI El El ■ In the RCC 
framework, the wave function of a single valence atom 
can be expressed as 



*<?>> = =e T {l + S n }\$ n ), 



(3.1) 



where $„) is the reference state constructed from the 
Dirac-Fock wave function |$ ) of the closed-shell config- 
uration by appending the corresponding valence electron 
as |<& n ) = a,Jj |$o) with a n representing addition of a va- 
lence electron n. Here T and S n are the RCC excitation 
operators which excite electrons from |$o) and \& n ), re- 
spectively. The amplitudes of these excitation are ob- 
tained by the following equations 



($ L \{H N e T }\* ) = 



($Z\{H N e T }S n \$ n ) 



+ (^\S n \<f n )AE n , 



(3.3) 



with the superscript L(= 1,2) representing the single and 
double excited states from the corresponding reference 
states and the wide-hat symbol over H^e T represent the 
linked terms of normal order atomic Hamiltonian Hn and 
RCC operator T. For the single and double excitations 
approximation (CCSD method), the corresponding RCC 
operators are denoted by 



T = T 1 +T 2 



and 



5i) 



Sin + S 2 



(3.4) 



(3.5) 



for the closed-shell and single valence configurations, re- 
spectively. Again, AE n in the above expressions is the 
corresponding valence electron affinity (negative of the 
ionization potential (IP)) energy which is evaluated by 



AE n = ($ n \{H N e T }{l + 5„}|$„). 



(3.6) 



In Eq. (|3.2p . we have considered only the single and 
double excitations, however we have incorporated con- 
tributions from important triple excitations (CCSD(T) 



method) perturbatively in Eq. (|3.3[) by defining 

■Sj^T* = HffT 2 -\- Hf^S 2n , (3-7) 
where the superscript pert denotes for the perturbation, 
and evaluating their contributions to AE n from these 
operators by 



AE^ = TlSlT. 



(3.8) 



After obtaining the amplitudes for T, the core excitation 
operator, we solve Eqs. (|3.3|) and (|3.6[) simultaneously to 
obtain the amplitudes for the S n operators. 

We extend the RCC ansatz for the perturbed atomic 
state in the presence of the electric dipole operator D by 
writing the total atomic wave function as 

|*n) = =e T+n {l + S n + k n }\^ n ), (3.9) 

where f2 and A„ are the first order corrections to the 
RCC operators T and S n , respectively. Since Eq. (|2.7|) 
is first order in the operator D, the above expression will 
reduce to 



l*n> 



e T {l + S n + fi(l + S n ) + A„}|$„}(3.10) 



(3.2) Now, separating the above wave function as Vt^) and 



I*! 1 )}, we get 



e T 



I n / 



e J Ml + S„)+A„}|$„>. (3.11) 



Following Eq. (|2.7p , we solve again the amplitudes for 
the modified operators as 



($ L \{H N e T Sl}\$ ) 
($„ L |{JJ^}A„|$„) 



-($ L |.De T |$ ) (3.12) 

-{^\{H^n(l + S n ) + De^ 
(l + 5„)}|$„) + ($^|A l ,|$„)A£; n , 

(3.13) 



where De T represents the connecting terms between D 
and T operators. Again in our CCSD approximation, we 
have 



and 



O = Oi + tt 2 



A n = A ln + A 2n . 



(3.14) 



(3.15) 



Therefore, the RCC expression for the dipole polariz- 
ability is given by 
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TABLE I: Comparison of dipole polarizabilities between different works in Sr + , Ba + and Ra + . 



— 

System 




nsi/2 


(n-l)d 3 / 2 (n-l)d 5 / 2 
a Qq al 


Mctnoas 


— — 

References 


Sr+(n = 


5) 














127.62 


145.86 -91.81 136.84 -116.02 


DF 


This work 






88.29(1.0) 


61.43(52) -35.42(25) 62.87(75) -48.83(30) 


CCSD(T) 


This work 






132.15 




HF 


[22] 






86.21 




non-rel. MBPT(2) 


[22] 






101.58 




non-rel. CCSD 


[22] 






97.91 




non-rel. CCSD(T) 


[22] 






121.33 




DK DF 


[22] 






79.89 




DK rel. MBPT(2) 


[22] 






94.31 




DK rel. CCSD 


[22] 






91.10 




DK rel. CCSD(T) 


[22] 






91.3(9) 


62.0(5) -47.7(3) 


LCCSD(T)+sum-over 


[16] 






89.88 


61.77 


Non-rel. +sum-over 


[23] 






93.3 


57.0 


Non-rel.+sum-over 


[24] 






84.6(3.6) 


48(12) 


Non-rel. +sum-over 


_ 






91.47 




Non-rel.+sum-over 


[25] 






86(11) 




Experiment + non-rel. 


[26] 


Ba + (n = 


= 6) 














184.49 


90.07 -45.07 87.66 -58.02 


DF 


This work 






124.26(1.0) 


48.81(46) -24.62(28) 50.67(58) -30.85(31) 


CCSD(T) 


This work 






213.47 




HF 


[22] 






110.60 




non-rel. MBPT(2) 


[22] 






148.24 




non-rel. CCSD 


[22] 






146.88 




non-rel. CCSD(T) 


[22] 






177.64 




DK DF 


[22] 






94.64 




DK rel. MBPT(2) 


[22] 






129.92 




DK rel. CCSD 


[22] 






123.07 




DK rel. CCSD(T) 


[22] 






124.15 




LCCSD(T)+sum-over 


[29] 






124.7 




Non-rel.+sum-over 


[25] 






126.2 




Non-rel.+sum-over 


[30] 






123.88(5) 




Experiment 


[27] 






125.5(10) 




Experiment 


[28] 


Ra+(n = 


= 7) 














164.66 


183.07 -114.70 143.77 -98.64 


DF 


This work 






104.54(1.5) 


83.71(77) -50.23(43) 82.38(70) -52.60(45) 


CCSD(T) 


This work 






257.00 




HF 


[22] 






123.23 




non-rel. MBPT(2) 


[22] 






186.23 




non-rel. CCSD 


[22] 






172.00 




non-rel. CCSD(T) 


[22] 






1/1^/17 




urv LT 


Tool 






79.80 




DK rel. MBPT(2) 


[22] 






110.48 




DK rel. CCSD 


[22] 






105.37 




DK rel. CCSD(T) 


[22] 






106.12 




CCSD(T)+sum-over 


[10] 






106.5 




LCCSD(T)+sum-over 


[32] 






106.22 




LCCSD(T)+sum-over 


[33] 


Abbreviations: HF — > Hartree-Fock. 
non-rel. — > non-relativistic. 
DK rel. —* scalar relativistic Douglas-Kroll method. 
MBPT(2) — » second order perturbation theory, 
sum-over — * sum-over intermediate states. 
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On 



<*1 0) |*1 0) > 

<$„|{1 + St }A{r>(l + S n ) + + (gnK4 +H+ ^)» f }A{l + S„}|$n) 

($„i{i + si}]vo{i + s„}i$„) 



(3.16) 



where we define L>i = (e T+ L>ie T ) and iVo = e Tt e T . The 

non-truncative series for Di and JVo are expanded us- 
ing the Wick's generalized theorem and truncated the 
series when the leading order non-accounted terms are 
below fifth order of Coulomb interaction. These oper- 
ators are then contracted with the fl to get fully con- 
tracted terms that give rise core electron contributions. 
The core-valence and valence correlation contributions 
are obtained from the open contraction between the op- 
erators with Q and f2{l+S„}+A n operators, respectively. 

Corrections due to the normalization of the wave func- 
tions are accounted by evaluating 



Norm = 



a radial grid given by 



n = r 



1 



(4.1) 



where i represents a grid point. The total number of grid 
points in our calculations is 750, the step size h is taken 
as 0.03 in the present case and tq is taken as 2 x 10~ 6 
atomic units. To construct the basis functions, we use 
Gaussian type orbitals (GTOs) defined as 



7 GTO 



in) 



(4.2) 



(^lAl^) + (^lAl^)] 1}, 

(3.17) 



where N n = ($„|{1 + Sl}N {l + S n }\$ n ). 



Here n K is the radial quantum number of the orbitals 
and on is a parameter whose value is chosen to obtain 
orbitals with proper behavior inside and outside the nu- 
cleus of an atomic system. Further, the a;,s satisfy the 
even tempering condition 



a f3 l 



(4.3) 



IV. RESULTS AND DISCUSSIONS 



A. General discussions 




(iv) H N 



FIG. 1: Break-down of RCC terms into lower-order pertur- 
bative diagrams. 



The orbitals used in the present work are generated on 



We have chosen same «o = 0.00525 and — 2.73 val- 
ues to construct the basis functions in Sr + , Ba + and Ra + , 
so that effects due to the different sizes of the systems 
can be compared. Finite size of the nucleus in these sys- 
tems are accounted by assuming a two-parameter Fermi- 
nuclear-charge distribution for evaluating the electron 
density over nucleus as given by 



p(n) = 



Po 



I -|_ e (n—c)/a ' 



(4.4) 



where c and a are the half-charge radius and skin thick- 
ness of the nucleus. These parameters are chosen as 



and 



a = 2.3/4(Zn3) 



5 2 7 

rms 



(4.5) 



(4.6) 



where r rms is the root mean square radius of the corre- 
sponding nuclei which is determined as discussed in [2l| . 

In Table H we present our DF and CCSD(T) results 
along with other available calculations and experimental 
results for the dipole polarizabilities of the ground and 
metastable states of Sr + , Ba + and Ra + . The differences 
between the DF and CCSD(T) results indicate the mag- 
nitudes of the electron correlation effects in the deter- 
mination of the dipole polarizabilities in these systems 
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TABLE II: Contributions from the DF and various CCSD(T) terms to the dipole polarizability calculations in Sr + . The 
subscripts c, n and v of the RCC terms correspond to the core, valence and virtual correlation contributions. D c and D n give 
the DF results from the core and valence orbitals. 



Terms 


c 2 ry 
OS 6i/ 2 

CX.Q 


4d D 3 / 2 


2 


Oil 


4d U5/2 

2 

QZq 




6.15 


6.15 


-0.25 


6.15 


-0.25 


£>„ 


121.47 


139.71 


-91.56 


130.69 


-115.78 


K 


4.98 


4.98 


-0.27 


4.98 


-0.27 


DM + cc 


0.10 


0.20 


-0.35 


0.41 


-0.41 


D v A ln + cc 


93.78 


67.05 


-38.70 


68.14 


-56.54 


D v A 2n + cc 


-2.87 


-2.54 


-0.86 


-2.43 


1.72 


Si v D v A ln + cc 


-3.74 


-6.06 


3.47 


-6.10 


5.05 


S2vD v A ln + cc 


-3.03 


-1.47 


0.77 


-1.39 


0.96 


S lv D v A2n + CC 


-0.04 


0.06 


0.03 


0.05 


-0.03 


S2vD v A2n + CC 


0.13 


-0.05 


0.03 


-0.06 


0.03 


Others 


0.13 


0.22 


-0.14 


0.24 


-0.15 


Norm 


-1.15 


-0.96 


0.60 


-0.97 


0.81 



using the CCSD(T) method. They are 45%, 48% and 
58% for the ground states of Sr+, Ba + and Ra + , respec- 
tively. They increase with the size of the system. How- 
ever, the correlation effects in the d 3 / 2 metastable states 
are 137%, 85% and 119% for the scalar polarizabilities 
and 159%, 83% and 128% for the tensor polarizabilities 
in Sr + , Ba + and Ra + , respectively. This shows that the 
correlation effects reduce in these states from Sr + to Ba + , 
but the presence of the core f-orbitals increases the corre- 
lation effects in Ra + . The correlation effects in the ^5/2 
metastable states are 117%, 73% and 75% for the scalar 
polarizabilities and 138%, 88% and 87% for the tensor 
polarizabilities in Sr + , Ba + and Ra + , respectively. This 
implies that the correlations in the d-metastable states 
do not depend upon the size but the internal structure of 
the systems. Our previous studies on the hyperfine struc- 
ture constants in these systems [l^EH, 15 1 had shown 
peculiar behavior of the core-polarization effects. These 
effects were comparatively smaller in the d 3 / 2 metastable 
states. In contrast, the correlation effects are larger in the 
d 3 / 2 metastable states compared to the d§/2 metastable 
states and the ground states in the dipole polarizabilities 
calculations. 

The upper limits to the error bars in these quantities 
were determined by taking the differences of the results 
obtained using CCSD(T) and CCSD methods and the 
inaccuracies due to the self-consistent results obtained at 
the DF levels by varying the number of GTOs considered 
in the calculations. These results are quoted inside the 
parentheses in Table fl] 

We explicitly present the diagrams in Fig. [T] corre- 
sponding to various RCC terms that are significant in 
determining the dipole polarizabilities. As seen from the 



figure, Fig. Ufa) which arises from the fully contracted 

terms of DCli corresponds to the core-correlation con- 
tributions. Its lower order terms corresponds mainly to 
the diagrams coming from the random-phase approxima- 
tion (RPA) . There are also core-correlation contributions 

arising from DU, 2l but they are relatively small and are 
not shown in Fig. [TJ The core-valence correlation con- 
tributions are determined by open diagrams from DQ,\ 
as shown in Fig. QJb). The most important correla- 
tion contributions arise through the valence correlation 
effects and they are shown in Fig. [ljc) and[TJd). Im- 
portant pair-correlation and core-polarization effects are 
accounted through DA ln , however core-polarization ef- 
fects arising from the perturbed doubly excited states 
are accounted through DA-2 n - The DF contributions in- 
volving the core, core- valence and virtual orbitals are the 

lowest order diagrams to the fully contracted Dfl\, open 
Dfti and -DA ln RCC terms, respectively. Based on the 
above mentioned correlation diagrams, we analyze their 
roles in different systems considered below. 



B. Sr+ 

There are no experimental results of the dipole polariz- 
abilities available for the ground and metastable excited 
states in Sr + . However, a number of calculations have 
been carried out using different methods and we have 
compared their results with the present work in Table 
HI Lim and Schwerdtfeger [22| have done comparative 
studies between the non-relativistic and scalar relativis- 
tic Douglas-Kroll calculations using four different many- 
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TABLE III: Contributions from the DF and various CCSD(T) terms to the dipole polarizability calculations in Ba + . The 
subscripts c, n and v of the RCC terms correspond to the core, valence and virtual correlation contributions. D c and D n give 
the DF results from the core and valence orbitals. 



Terms 


a~ 2 o 
OS 6i/ 2 


Oil 


5d D 3/ 2 

2 


i 

(Xq 


5d D 5 / 2 

2 




11.73 


11.73 


-0.46 


11.73 


-0.46 


£>„ 


172.76 


78.33 


-44.61 


75.93 


-57.56 


K 


9.35 


9.35 


-0.56 


9.35 


-0.56 


DM + cc 


0.23 


0.33 


-0.64 


0.82 


-0.82 


D v A ln + cc 


133.01 


49.20 


-25.61 


50.17 


-36.23 


D v A 2n + cc 


-4.93 


-3.23 


-1.36 


-3.05 


2.19 


SivD v A ln + cc 


-6.58 


-4.45 


2.39 


-4.42 


3.37 


S2vD v Ai„ + cc 


-5.18 


-1.98 


0.84 


-1.83 


0.89 


S lv D v A2n + CC 


-0.06 


0.05 


0.03 


0.03 


-0.01 


S2vD v A2n + CC 


0.27 


-0.10 


0.05 


-0.09 


0.06 


Others 


0.21 


0.45 


-0.25 


0.49 


-0.33 


Norm 


-2.06 


-0.81 


0.49 


-0.80 


0.59 



body methods. They demonstrate the importance of the 
relativistic methods to calculate dipole polarizabilitics. 
Jiang et al. [l6| have used El matrix elements obtained 
using the linearized RCC method with the singles, dou- 
bles and partial triple excitations (LCCSD(T)) to eval- 
uate the valence correlation contributions for a few in- 
termediate states. The core-correlations are accounted 
through the RPA method and contributions from higher 
states were estimated using the DF method. Mitroy et 
al. [HI have used a non-relativistic method using the 
sum-over-states approach to determine polarizabilities of 
the ground and d-state. As seen in Table UJ the dipole 
polarizabilities of the 4d3/2 and 4ds/2 states are not the 
same and they cannot be evaluated separately using a 
non-relativistic method. However, our ground state po- 
larizability for Sr + agrees with their result. Similar ap- 
proaches were also employed by Barklem and OMara [2J] . 
Patil and Tang [25| have employed a summation and in- 
tegration approach to determine the ground state polar- 
izability. Recently, Nunkaew et al. [2f| have estimated 
El matrix elements using the non-relativistic theory and 
microwave resonance measurements in Sr and have ex- 
tracted dipole polarizability of the ground state of Sr + . 

In Table [Til we present the individual contributions 
from RCC terms to the dipole polarizability calculations 
in Sr + . Our core-correlation contributions are 4.98 au 
and —0.27 au for the scalar and tensor dipole polariz- 
abilities, respectively. Clearly, the CCSD(T) result for 
the scalar dipole polarizability is smaller than the pre- 
viously estimated values. On the otherhand, the core- 
correlation to the tensor polarizability vanishes in the 
non-relativistic theory, but it is finite in our approach, 



although small in magnitude. Jiang et al. [lq have ne- 
glected this contribution in their calculations. We have 
also given DF results from the core (Z? c ) and virtual 
(D v ) orbitals separately in the same table. Our DF re- 
sult and that reported by Lim and Schwerdtfeger [22j 
differ. Comparing our DF results given in Table UJ and 
Table ILT1 it seems that Lim and Schwerdtfeger have not 
included core correlation contributions at the DF level. 
Again, the lowest order contributions to D c and D v Ai n 
terms correspond to D c and D Vl respectively. The dif- 
ferences between the lowest order and all order results 
seem to be significant in this system. The largest contri- 
butions to the final results come from D„Ai„ as it con- 
tains DF results due to virtual orbitals in it. Contri- 
butions from D v A2 n correspond to doubly excited per- 
turbed states and they are also large in both the ground 
and metastable states. Therefore, the exclusion of these 
contributions in the sum-over-state approach may not be 
appropriate. Again, normalization corrections (Norm) 
are also non-negligible. 



C. Ba+ 

Two experimental results with small uncertainties 
[27L |28| are available for the ground state dipole po- 
larizability in Ba + . There have also been studies of 
this quantity by Lim and Schwerdtfeger [22] ], Iskrenova- 
Tchoukova and Safronova [29[ have employed El matrix 
elements from the LCCSD(T) method in the sum-over- 
states approach using a few states for the valence correla- 
tion effects and estimating the core-correlation and core- 
valence correlation contributions from lower order pertur- 
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TABLE IV: Contributions from the DF and various CCSD(T) terms to the dipole polarizability calculations in Ra + . The 
subscripts c, n and v of the RCC terms correspond to the core, valence and virtual correlation contributions. D c and D n give 
the DF results from the core and valence orbitals. 



Terms 


7s 6i/ 2 

CX.Q 


ok 


6d D 3/2 

2 


! 


6d D 5/2 

2 




15.56 


15.56 


-0.56 


15.56 


-0.56 


£>„ 


149.10 


167.51 


-114.14 


128.21 


-98.07 


K 


11.66 


11.66 


-0.71 


11.66 


-0.71 


DM + cc 


0.60 


0.21 


-0.54 


1.03 


-1.03 


D v A ln + cc 


107.74 


91.30 


-54.32 


85.59 


-62.17 


D v A 2n + cc 


-4.15 


-5.92 


-2.84 


-4.85 


3.68 


Si v D v A ln + cc 


-5.17 


-7.62 


4.51 


-6.65 


5.13 


S2vD v A ln + cc 


-4.93 


-4.64 


2.68 


-3.44 


1.76 


S\ v D v A.2n + CC 


-0.02 


0.05 


0.05 


0.02 


-0.01 


S2vD v A2n + CC 


0.25 


-0.18 


0.07 


-0.15 


0.10 


Others 


0.29 


0.61 


-0.34 


0.64 


-0.43 


Norm 


-1.73 


-1.76 


1.21 


-1.47 


1.08 



bation theory to determine this quantity. Other available 
calculations (2f| [3(| are based on non-relativistic meth- 
ods. Again, there are no other results available for the 
metastable d-states in Ba + to compare with our results. 
However, we have also carried-out a sum-over-states cal- 
culation using the El matrix elements from the CCSD(T) 
method [3l| that agrees with our ab initio results. 

We present contributions from individual RCC terms 
to the dipole polarizabilities calculations on Ba + in Table 
IIIII The trends of these correlation effects seem to be the 
same as in Sr + . However, the core correlation effects in 
this system seem to be almost twice than in the case of 
Sr + . The core- valence correlations coming through the 
open D v fl diagrams are also larger than Sr + . Contribu- 
tions from the doubly excited perturbed states and cor- 
rections due to the normalization of the wave functions 
also seem to be significant. 



d-metastable states using the sum-over-states approach 
with the El matrix elements from CCSD(T) method and 
approximated core-correlation and core-valence correla- 
tion effects. 

In Table HVl we present contributions from individual 
RCC terms to these results. The trend of the correla- 
tion effects in the ground state seems similar to those of 
Sr + and Ba + , but due to the presence of core f-electrons, 
the behavior of the correlation effects is a little different 
for the metastable d-states. The size of core-correlation 
is slightly larger than that of Ba + , but the difference is 
not as large as it was between Sr + and Ba + . In contrast 
to Ba + where the ab initio and sum-over-states results 
match, we found discrepancies in this system. The dis- 
crepancies are mainly because of the inclusion of the dou- 
bly excited states in the present work, but there could be 
cancellations in Ba + due to which the discrepancies are 
small. 



D. Ra+ 

There are also no experimental results available for the 
dipole polarizabilities in Ra + . In the same work as men- 
tioned above, Lim and Schwerdtfeger [12] have also cal- 
culated this quantity in the ground state of Ra + u sing 
various many-body methods. Safronova et al (32l . 1 3 3l ] 
have also evaluated this result using the sum-over-states 
approach. Their valence correlation effects are evalu- 
ated using El matrix elements for a few important states 
from the LCCSD(T) method and core-correlation and 
core- valence correlations are evaluated using lower order 
many-body methods. In our earlier work [lfj, we had 
also evaluated dipole polarizabilities in the ground and 



E. Applications to the optical clocks 

All the ions considered in this work are important can- 
didates for optical clocks 0, M, B B E3| • There has been 
an absolute frequency measurement of the 5s 2 5i/2 — > 
4d 2 D 5 / 2 transition in 88 Sr+ by Madej et al. [5]. One 
of the largest uncertainties due to the applied electric 
field comes from the quadratic Stark shift. In fact, this 
shift was earlier over estimated due to the large error 
bars in the calculated dipole polarizabilities of the 5s 
2 Si/2 and 4d 2 D 5 / 2 states. Madej et al. had used 
a i(5s 1/2 ) = (1.40±0.06) x 1CT 39 C 2 s 2 kg" 1 where we ob- 
tain this result as (1.46 ±0.02) x 10~ 39 C 2 s 2 kg" 1 . The 
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scalar and tensor polarizabilities of the 4d 2 -D 5 / 2 were 
used in [| as aJ(4d 5/2 ) = (8 ±2) x lO" 40 C 2 s 2 kg" 1 and 
a 2 (4d 5/2 ) = (-7±2) xlO~ 40 C 2 s 2 kg" 1 , respectively. We 
obtain these results as aJ(4d 5/2 ) = (10.37±0.12) x 10^ 40 
C 2 s 2 kg- 1 and a§(4d 5/2 ) = (-8.05 ± 0.05) x lO" 40 
C 2 s 2 kg -1 . Using Eq. (|2.9j) and our results, we ob- 
tain the shift rate, which is defined as 7 = ^^/^ , of 
the 5s 2 5 1/2 state as (1.10 ± 0.01) ^Hz/(V/m) 2 against 
(1.06 ± 0.04) ^Hz/(V/m) 2 of Madej et al. Similarly, we 
obtain 7 = (-0.78 ±0.02) ^Hz/(V/m) 2 against results of 
Madej et al. as 7 = (-0.6 ± 0.2) /iHz/(V/m) 2 in the 4d 
2 D 5 / 2 state using only the scalar polarizability. However 
assuming the direction of the electric field lies in the z- 
direction, we obtain 7 = (-1.27 ± 0.03) /iHz/(V/m) 2 , 
7 = (-0.91 ± 0.02) /iHz/(V/m) 2 and 7 = (-0.18 ± 
0.01) MHz/(V/m) 2 for M = 1/2, M = 3/2 and M = 5/2, 
respectively. 

Using Eq. (|2. 14[) and the above results, we also obtain 
the black-body radiation shift at T = 300K in the 5s 
2 Si/ 2 -> 4d 2 L> 5/2 transition in 88 Sr+ as (0.22 ± 0.01) 
Hz and that is an improvement of 10% over the result of 
Madej et al. @]. 

It appears that both 137 Ba + and 138 Ba + will be suit- 
able candidates for an optical clock || [1J , but each has 
some advantages and disadvantages in controlling the 
systematic errors. For the 6s 2 Si/ 2 — > 5d 2 D 5 / 2 tran- 
sition in 138 Ba + , it would be possible to use techniques 
similar to the measurement of the frequency in the optical 
transition in 88 Sr + mentioned earlier. However, one has 
to encounter the electric quadrupole shift in the 5d 2 D 5 / 2 
state for this case. It is possible to overcome this particu- 
lar shift by considering the possible F — 2(6si/ 2 ) — ► F = 
0(^3/2) hyperfine transition in 137 Ba + . In this transi- 
tion, one has to again estimate the possible quadratic 
Zeeman shifts because of finite nuclear magnetic and 
quadrupole moments. Our dipole polarizability for the 
6s 2 S 1/2 state is given by (2.05 ± 0.02) x 10~ 39 C 2 s 2 
kg -1 . The scalar and tensor polarizabilities of the 5d 
2 L» 3/2 state are given by (8.05 ±0.07) x 10~ 40 C 2 s 2 kg" 1 
and (-4.06±0.05) x 10" 40 C 2 s 2 kg" 1 , respectively. Simi- 
larly, the scalar and tensor polarizabilities of the 5d 2 D 5 / 2 
state are given by (8.35 ± 0.09) x 10~ 40 C 2 s 2 kg" 1 and 
(-5.09 ± 0.05) x 10~ 40 C 2 s 2 kg" 1 , respectively. Due to 
the choice of the hyperfine transition in 137 Ba + , the ten- 
sor polarizabilities of these states are zero and hence the 
polarizabilities of the atomic and hyperfine states are the 
same. The shift rates are (-1.55 ± 0.01) ^iHz/(V/m) 2 
and (-0.61 ± 0.01) /iHz/(V/m) 2 in the 6s 2 S' 1/2 and 
5d 2 D 3 / 2 states, respectively. For 138 Ba + , by consid- 
ering particular M values of the 5d 2 D 5 / 2 state and as- 
suming that the electric field lies in the z-direction, we 
can evaluate the Stark shifts. They are 7 = (—0.94 ± 
0.02) ^Hz/(V/m) 2 , 7 = (-0.71 ±0.02) /m Z /(V/m) 2 and 
7 = (0.14 ± 0.01) /iHz/(V/m) 2 for M = 1/2, M = 3/2 
and M — 5/2, respectively. As can be noticed, the result 
for M = 5/2 has opposite sign than other M values. The 



Stark shifts in these states can be easily estimated using 
these results for a given applied electric field. 

The black-body radiation shift at T = 300K in the 6s 
2 S1/2 — > 5d 2 Dc ) / 2 transition in this system is given as 
(0.64 ±0.12) Hz. 

Similarly as we had reported earlier [lOj, both 223 Ra + 
and 225 Ra + have the same advantages like 137 Ba + and 
138 Ba + , respectively, for considering as optical clock can- 
didates. In fact, all the low-lying energy levels in these 
ions are in optical region which will be an advantage for 
the experimentalists to measure the 7s 2 5' 1 / 2 — > 6d 2 D 3 / 2 
or 7s 2 S 1/2 -> 6d 2 D 5/2 or F — 2(7s 1/2 ) -» F = 0(6d 3/2 ) 
transition frequencies more precisely than other candi- 
dates. Recently, 213 Ra whose half-lifetime is around 2.75 
m was produced at KVI [34| in the accelerator method 
and its single ion shares the same advantage with 225 Ra + 
for becoming suitable candidate for the optical clock. 
Now assuming that due to the suitable choice of hyperfine 
states in 223 Ra + [10] like the case for 137 Ba + , the tensor 
polarizability contribution to the Stark-shift will be zero 
and hence using our dipole polarizability results, we ob- 
tain the Stark shift rates as (-1.31 ± 0.02) /xHz/(V/m) 2 
and (-1.05 ±0.02) ^Hz/(V/m) 2 in the 7s 2 S' 1/2 and 6d 
2 D 3 /2 states, respectively. For other isotopes discussed 
above, by considering particular M values of the 6d 2 D 5 / 2 
state and assuming that the electric field lies in the z- 
direction, the Stark shifts are evaluated as 7 = (— 1.56± 
0.03) /xHz/(V/m) 2 , 7 = (-1.16 ±0.02) f iIlz/(V/m) 2 and 
7 = (0.29 ± 0.01) /iHz/(V/m) 2 for M = 1/2, M = 3/2 
and M = 5/2, respectively. The result for M = 5/2 has 
opposite sign than other M values like in 138 Ba + . There- 
fore, the Stark shifts in Ra + can be estimated accurately 
using our results for a given applied electric field. 

The black-body radiation shift at T = 300K in the 7s 
2 Si/2 — > 6d 2 F)^j2 transition in Ra + is given as (0.19 ± 
0.02) Hz. 

From the above Stark shift ratios and BBR shifts in 
the considered ions, it is found that these systematic er- 
rors are small in Ra + which further supports along with 
its energy level locations that it will be one of the most 
suitable candidates for optical clock. In fact, a possi- 
ble atomic clock with uncertainty in the order of 10~ 
seems feasible from these results along with the prelimi- 
nary analysis of Doppler's shifts [35j in Ra + . 

V. CONCLUSION 

We have employed the relativistic coupled-cluster 
method to determine ab initio results for the dipole po- 
larizabilities of the ground and the metastable d-states in 
the singly ionized strontium, barium and radium. Elec- 
tron correlation effects arising through various coupled- 
cluster terms are given individually and comparative 
studies are performed for these three ions. Using the 
results we have obtained, Stark shifts and black-body 
radiation shifts for these ions are estimated. Using our 
results, we were able to reduce the errors of the measured 
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frequency for the optical clock in 88 Sr + . Our calculations 
of the Stark and black-body radiations shifts in both Ba + 
and Ra + could be used to remove the systematic errors 
in the proposed optical clock experiments for these ions. 
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